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OBJECTIVE 


Develop  an  algorithm  and  subroutine  for  evaluating  J  (z)  to  be  used  in 

N 

determining  the  scattering  of  plane  elastic  waves  from  a  single,  hollow  cylin¬ 
drical  cavity.  Optimal  accuracy  and  computation  time  are  both  essential. 

RESULTS 


1.  Two  series  were  developed  for  computational  purposes:  one  for  the 

real  part  of  J  (z)  and  one  for  the  imaginary  part  of  J  (z). 

N  N 

2.  A  subroutine  was  developed  for  evaluating  J  (z)  by  utilizing  double- 

N 

precision  arithmetic.  Included  in  the  output  is  the  complex  number  EZ  where 
the  real  and  imaginary  parts  represent  the  magnitudes  of  upper  bounds  of 
errors  in  respective  parts  of  J„(z). 

3.  Computational  studies  strongly  suggest  that  for  |z|  £  10,  the  magni¬ 

tudes  of  the  errors  in  the  real  and  imaginary  parts  of  JM(z)  are  less  than 
10  .  From  the  same  studies,  computational  times  were  always  less  than  0.04 

second . 


RECOMMENDATIONS 

1 .  Although  the  motivation  for  designing  the  algorithm  and  subroutine  for 
computing  J^Cz)  was  for  use  in  the  specified  scattering  problem,  satisfactory 
usage  in  a  wide  range  of  scientific  problems  is  projected. 

2.  Subroutine  JNZ  can  also  be  used  for  Bessel  functions  of  real  arguments 
by  setting  the  imaginary  part  of  z  to  zero.  Computational  studies  show  a  high 
degree  of  accuracy  (magnitude  of  the  maximum  error  less  than  10-6)  for  z 
taking  on  integer  values  from  1  to  10. 
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INTRODUCTION 


Computational  methods  for  evaluating  Bessel  functions  of  the  first  kind 
have  been  a  necessary  part  of  scientific  research  and  technology  for  gener¬ 
ations.  Exhaustive  tables  [1]  have  been  published  for  real  arguments,  while  a 
similar  publication  for  complex  arguments  [2]  is  much  more  limited  in  scope. 

A  recent  need  arose  at  NOSC  for  these  functions  with  complex  arguments  in 
determining  the  scattering  of  plane  elastic  waves  from  a  single  hollow  cylin¬ 
drical  cavity.  To  meet  this  need,  priorities  of  accuracy  and  computational 
time  were  set  as  number  one  and  number  two  goals,  respectively. 


ALGORITHM 


Although  some  analysts  will  question  the  efficiency  of  the  series  repre¬ 
sentation  for  computational  purposes  (both  in  terms  of  accuracy  and  computa¬ 
tional  time),  this  study  shows  that  the  approach  does  reasonably  well.  Of 
course,  no  universal  algorithm  has  been  designed  to  successfully  compute  J„(z) 
for  the  entire  domain  of  existence  [3]  of  the  complex  numbers. 

The  series  representation  of  the  Bessel  function  is  given  by: 

»  12 
,  N  ^  (-  ±2*) 

JN<Z)  =  (2Z)  2^  Kl  (N+K)  I 

K=0 

=  U  +  iV, 

for  z  =  p  (cos  0  +  isin  0). 


1.  Harvard  Computation  Laboratory  Annals,  vol  III-IX,  Harvard  University 
Press,  Cambridge,  MA,  1947-48. 

2.  Tables  of  the  Bessel  Functions  Jg(z)  and  J^z)  for  Complex  Arguments, 
Columbia  University  Press,  New  York,  1943. 

3.  Private  consultation  with  Dr  Gordon  Martin. 
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K=0 


2  r  (K+i )  r  (n+k+i) 


V 


z 


K  2K 

(-1)  p  sin  (Nt2K)  9 
2K 

2  T  (K+1)  T  (N+K+1) 


Thus  we  have  two  independent  series  representing  the  real  and  imaginary 
parts  of  the  Bessel  function. 

COMPUTATIONAL  RESULTS  AND  CONCLUSIONS 

All  programming  for  this  report  was  done  in  ASCII  Fortran  for  the  UNIVAC 
1100/82. 


The  subroutine  for  computing  JN(z)  is  referenced  by  the  statement 
CALL  JNZ  ( JN , NP 1 ,  Z,  EZ) 

where 

will  contain  the  double-precision  complex  output  JN(z) 

NP1  is  the  integer  order  of  the  Bessel  function 
Z  is  the  double-precision  complex  argument 

EZ  is  a  double-precision  complex  number  with  the  real  and  imaginary  parts 

being  upper  bounds  for  the  magnitudes  of  respective  parts  of  the  relative 

error  in  J  . 

N 

Subroutine  JNZ  carries  out  the  computations  in  double-precision  arithmetic 
and  uses  double-precision  system  subroutines: 

DREAL,  DIMAG,  DSQRT,  DATAN2 ,  DSIN,  DCOS,  DGAMMA,  DCMPLX. 


Since  reference  1  was  available  for  comparison  of  computations  with  real 


arguments,  it  was  decided  to  compute  JN(z)  for  N  =  1 ,  2,  . ..,  27,  and  z  =  1, 

2,  3,  ...,  15.  In  general,  the  comparisons  showed  the  magnitudes  of  the 

—8 

differences  to  be  less  than  10  in  the  worse  cases  [12  values  of  J^Cz)]  for 

z  =  1,  2,  3,  4,  and  5.  For  N  >_  9,  the  computations  agreed  for  the  10  places 

of  output  with  reference  1.  For  z  =  6,  7,8,9,  and  10,  two  values  showed 

—8  —7 

comparison  differences  greater  than  10  but  less  than  10  .  All  other  dif- 

«8 

ferences  were  significantly  less  than  10  .  For  these  same  cases  (z  =  6  to 

10),  the  computations  agreed  with  reference  1  for  the  10  places  of  output  for 
N  14.  The  computations  for  z  *  11,  12,  13,  14,  and  15  showed  significantly 
larger  differences  than  the  other  computations. 

In  general,  whenever  differences  were  noted  between  reference  1  and  the 
computations  for  z  -  1  to  10,  the  magnitudes  of  the  computations  were  less 
than  the  corresponding  magnitudes  of  reference  1.  For  z=  11  to  15,  the 
computational  values  were  algebraically  less  than  corresponding  values  given 
in  reference  1 . 


Tables  1  through  26  represent  some  results  for  selected  values  for  the 
argument  z,  where  'N'  is  the  order  of  JN»  JN  is  the  computed  complex  values  of 
the  Bessel  function,  and  the  time  of  computation  (TIME)  is  given  in  seconds. 
Listed  under  the  heading  BESL.EQ.,  are  the  complex  values  of  the  Bessel  equ¬ 
ation  for  the  specified  as  its  solution.  (Note  that  the  value  of  the 
Bessel  equation  should  be  zero  when  satisfied. )  Listed  under  the  heading, 

REL. ERROR,  is  the  output  '  EZ'  from  the  subroutine  JNZ. 


The  following  expression  is  the  Bessel  equation  given  in  terms  of  its 


solution  J  (z). 
N 


2  r 

f  [v 


2'*>  -*V*>  -fw* 


*  W2’] 


+  (z2-N2)  J  (z)  =  0. 


Therefore,  the  values  listed  under  the  heading,  BESL.EQ.,  represent  'com¬ 
plex  zeros'  and  along  with  values  listed  under  REL. ERROR  give  some  indications 
of  the  quality  of  the  computations.  From  this  standpoint,  tables  16,  17,  and 
18  are  the  worst  results. 

Although  the  number  of  computations  carried  out  here  has  not  been  exhaus¬ 
tive,  it  is  conceivable  that  for  |z|  £  10,  JNZ  can  carry  sufficient  accuracy 
for  many  needs.  Additional  computational  studies  strongly  suggest  that  for 

| z |  £  yVF,  the  magnitudes  of  the  real  and  imaginary  parts  of  both  EZ  and  the 

-9 

error  in  the  computational  value  of  the  Bessel  equation  are  less  than  10  . 
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